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ABSTRACT 



Context. The light echo phenomenon that accompanied the 2002 eruption of V838 Mon allows one to study the properties of the 
diffuse dusty matter in the vicinity of the object. 

Aims. We are aiming at obtaining estimates of the optical thickness of the circumstellar matter in front of V838 Mon, as well as 
optical properties of dust grains in the echoing medium. In particular, we are interested in studying whether the echoing medium can 
be responsible for the observed faintness of the B-type companion of V838 Mon when compared to three B-type stars that are seen in 
the vicinty of V838 Mon and are believed to be at the same distance as V838 Mon. 

Methods. We used the V838 Mon light echo images obtained by the Hubble Space Telescope (HST) in different filters and epochs. 
From the images we derived the total brightness of the echo and its surface brightness. The results of the measurements were compared 
to model light echoes. 

Results. The present study allowed us to estimate the optical thickness of the matter in front of the object and the mean cosine value 
of the scattering angle of dust grains in three HST filters. 

Conclusions. The optical thickness of the echoing matter is not sufficient to explain the observed difference in brightness between the 
B-type companion of V838 Mon and the other three B-type stars observed in the vicinity of V838 Mon. Implications of this result are 
discussed. Our estimate of the mass of the diffuse matter seen in the light echo shows that the matter cannot have resulted form a past 
mass loss activity of V838 Mon. We probably observe remnants of an interstellar cloud from which V838 Mon and other members of 
the observed cluster were formed. 

Key words, stars: individual: V838 Mon - stars: peculiar - stars: late-type - circumstellar matter 
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1. Introduction 

V838 Monocerotis was discovered as an eruptive star at the 
beginning of January 2002. The main eruption started at the 
beginning of F ebruary 2002 however, and lasted about two 
months (see e.g. Munari et"al]l2002al : iKimeswenger et al]|2002t 



ICrause etal ] I2003h . The main characteristic of the eruption, 
which is different in V838 Mon from other stellar eruptions (e.g. 
classical novae), was that the object evolved to progressively 
lower effective temperatures a nd declined as a very late M-type 
supergiant (e. g. lTvlendall2005h . 

ITvlenda & Sokeii ( |2006|) showed that the eruption of V838 
Mon cannot be explained by a thermonuclear runaway similar 
to classic al novae, nor by a late H e-shell flash. They proposed, 
following [Soker & Tvlendal (l2003h . that the event resulted from 
a merger of two stars. This idea obtained strong support from 
an analysis of archive observations of the progenitor of VI 309 
Sco. The eruption of this object was obser ved in 2008 and wa s 
of the same t ype as that of V838 Mon dMason et all l2010h . 
ITvlenda et al] (1201 lal) showed that the event resulted from a 
merger of a contact binary. 



Send offprint requests to: R. Tylenda 

* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
Science Institute. STScI is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract NAS 5-26555. 



The eruption of V838 Mon was accompanied by a spectac- 
ular light echo. A series of fantastic images of the light echo 
have been obtained by Hubble Space Telescope (HST, see e.g. 
IBond et al. 2003;. Bond 2007) . Analyses of the light echo's struc- 
ture and properties allowed ast rophysicists to study the dusty 
medium around V838 Mon (e.g.'Tvlenda'2004^ and to determine 
the distance to the object (Sparks et al. 2008). 

V838 Mon has a spe ctroscopic companion, a B-type star dis- 
covered bv iMunari et a l. (2002b) in September-October 2002, 
when the eruption remnant became so cool that it almost disap- 
peared from th e optical. The spectru m of the companion, classi- 
fied as B3 V bv lMunari et al ] (l2002bh . was seen undisturbed until 
the fall of 2006, when an eclipse-like event was observed (e.g. 
Munari et al. 2007a). The companion was eclipsed for about 
two months. In 2004/2005, emission lines of [Fell] appeared 
in the spectrum of V838 Mon. The lines strengthened with time 
and reached a max imum at the time of the eclipse-like event. 
ITvlenda et al.l (120091) showed that the profiles and evolution of 
the [Fe II] emission lines, as well as the eclipse-like event re- 
sulted from interactions of the B-type star radiation with the mat- 
ter ejected by V838 Mon in the 2002 eruption. After reappear- 
ance from the eclipse-like event, the companion began to fade 
again and in 2008 traces of the B-type spectrum, as well as of the 
[Fell] emission lines, disappeared. Apparently the companion 
became completely em bedded in the dusty ejecta of V838 Mon 
(ITvlenda etalJl2011bh . 
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lAfsar & Bondl (l2007h found a sparse, young cluster, contain- 
ing three B-type main-sequence stars, in the field of V838 Mon. 
The cluster is at a distance of ~6.2 kpc, wh ich is very close to the 
value of -6.1 kpc derived bv lSparks et all (2008) for V838 Mon 
from their analysis of the polarization of the light echo. The red- 
dening of the cluster and that of V838 Mon and its companion 
are also very similar {Eb-v - 0.84 versus 0.9, respectively). It 
seems very likely that V838 Mon and its B-type companion are 
also members of the cluster. 

However, there is a problem with the observed brightness of 
the B-type companion of V838 Mon wh en compared to th e mag- 
nitudes of the B-type stars of the cluster. iMunarieTaT (2005) de- 
rived V - 16.05 + 0.05 for the companion of V838 Mon. Using 
the compilation of photometric data of V. GoranskiJ^, one de- 
rives V =16.19 + 0.03 from the data obtained in September- 
November 2002. The B-type companion is therefore fainter by 
1.2-1.4 mag when compared to s t ar 9 (classified as B3 V) in 
the cluster of lAfsar & Bond (l2007h . lAfsar & Bondl (l2007h inter- 
preted this as evidence that the companion is partly submerged 
in the dusty ejecta of V838 Mon. However, if a dusty medium 
results from a stellar eruption, it is very unlikely that it is ho- 
mogeneous. Therefore, if the ejecta flows in front of a star, one 
expects to observe a high variability of the stellar light. This was 
not observed in the case of the B-type companion of V838 Mon 
until the fall of 2006. As can be seen from the light curve of 
V838 Mon between fall 2002 and fall 2006 (see e.g. the data of 
Goranskij'), the U brightness (dominated by the B-type com- 
panion in this epoch) remained constant within ±0.1 mag. The 
light curve exhibit s the same behavio ur in the B band after the 
2002 eruption (see Munari et al. 2005). As discussed above, in- 
teractions of the V838 Mon ejecta with the companion started at 
the earliest in the fall of 2006 and then indeed resulted in a high 
variability of the observed light of the companion, eventually 
leading to the complete disappearance of the star 

The problem of the observed brightness of the B-type com- 
panion remains unresolved. A possible explanation is that V838 
Mon and its companion suffer from additional, local extinction 
of AAv - 1.3, compared to the three B-type stars of the clus- 
ter. A standard extinction of this value would imply additional 
reddening of AEb-v - 0.4. This is significantly higher than the 
observed difference between the reddening of the V838 Mon and 
its companion, and that of the cluster (AEb-v $0.1, see above). 
This implies that the local extinction, if present, would have to 
be almost grey in the optical, i.e. R - Av/Eb-v ~ 10, which 
would additionally imply local dust dominated by large grains 
(a k 0.3 yum). 

This idea seems to be conceivable. Young stellar clusters are 
usually known to contain a significant amount of diffuse matter 
consisting of remnants of a more massive complex, from which 
the cluster was formed. For V838 Mon there is observational 
evidence that diffuse matter is indeed present in the near vicin- 
ity of the object. One piece of evidence is the light echo. Radio 
observations in the CO rotational lines also show the presence 
of molecular matter, m ost probably related to t he du sty medium 
seen in the light echo dKamihski et alj|2007l 1201 lb . The fight 
echo offers possibilities to study parameters and properties of 
the scattering matter. This is the aim of the present study. 

As we show in detail below, an analysis of the brightness 
of the echo at different wavelengths and its evolution with time 
allows us to estimate several parameters of the scattering mat- 
ter. This includes the optical thickness of the matter, as well as 
dependeces of the scattering coefficient on the wavelength and 



' http://jet.sao.ru/jet/~goray/v838mon.htm 



the scattering angle. In this way we can verify whether the lo- 
cal extinction can be responsible for the relative faintness of the 
B-type companion of V838 Mon when compared to the other 
B-type stars of the cluster The results of the echo study also al- 
low us to derive constraints on the nature of the dust grains and 
on the mass of the scattering material. The latter is important for 
discussing the nature of the matter, i.e. whether it is of interstellar 
origin or rat her, a s suggested by some authors (e.g. Bond et al] 
2003; B ondll2067h . that it resulted from past mass loss activity 
of the progenitor of V838 Mon. This can have consequences for 
the nature of the progenitor itself and the mechanism of the 2002 
eruption. 

2. Observational material and its reduction 

We have used the archival images of V838 Mon's light echo ob- 
tained with HST in 2002 on April 30, May 20, September 2, 
October 28, and December 17 (HST proposals 9587, 9588, and 
9694). The first image was taken in the F435W filter only. The 
other ones were obtained in the F435W, F606W, and F814W 
filters. The observations were performed with the Advanced 
Camera for Surveys (ACS) combined with the detector of the 
Wide Field Channel (WCS). Some of the images were taken in 
the polarimetric mode, the rest in the standard (non-polarimetric) 
mo de. The o bservations have been described in more detail in 
Sparks et al.l (12008) . The pipeline-reduced data, which we ex- 
tracted from the HST archive, were further processed using the 
multidrizzl^ package (version 3.3.8). The non-polarimetric ob- 
servations for a given date and filter were combined in multidriz- 
zle to produce images in units of counts per second. The po- 
larimetric observations were obtained with three retarder angles 
(0°, 60°, 120°). We first combined frames with the same retarder 
angle for a given filter and date. Then the polarimetric images 
were scaled using calibration corrections taken from ACS Data 
Handbook Finally the data were combined into the Stokes / 
(total intensity) image in each filter 

The following reduction was performed using StarHnliQ 
packages, mainly with KAPPA and GAIA. We measured the 
background sky level and its variations in each image by per- 
forming statistics on a large number of pixels free of stellar and 
echo emission and instrumental defects. After correcting images 
for the sky level, those measurements were used to blank all pix- 
els below 1.5cr of sky variations and above the maximum echo 
brightness for a given date. From the resulting images, we man- 
ually removed the remaining parts of stellar diffraction patterns 
and cosmic -ray hits. Finally, we measured the light echo region. 

3. Measurements 

The results of the measurements are given in Table [T| The up- 
per, middle, and bottom parts of the table present the results 
obtained with the F435W, F606W, and F814W filters, respec- 
tively. The dates of the observations are given in column (1), 
while column (2) notes the observing mode used to obtain the 
images. Total count rates from the echo are given in column (3). 
Column (4) gives the total number of pixels, which registered a 
measurable signal from the echo. A total flux from the light echo 
can be calculated by multiplying the total count rate number by 
the photflam conversion factor (Sirianni et al.l |20()5 ). which is 
equal to 3 .08 x 1 0" ' 7 .725 x 1 0"'^ and 6.94 x 1 0" W m"- fim- ' 

' http://stsdas.stsci.edu/multidrizzle 

^ http://www.stsci.edu/hst/acs/documents/handbooks/currentDHB 
* http://starlink.jach.hawaii.edu/starlink 
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per count s ' for the F435W, F606W, and F814W filters, respec- 
tively. The resulted fluxes are presented in column (5). The ACS 
pixel scale is (X/OS x (X/OS, which allows one to easily calculate a 
mean surface brightness, ^b, from the mean count rate summed 
over 400 pixels (corresponding to 1 arcsec^) and the photflam. 
The results are listed in column (6). 

The data in columns (3)-(6) refer only to the part of the 
light echo surface that was not contaminated by field stars and 
V838 Mon itself. The procedure of removing diffraction patterns 
of field stars introduced uncertainties in the measurements. This 
is particularly the case for the total flux, which is expected to 
be systematically underestimated in column (5) of Table [1] The 
effect is difficult to estimate, as each time different parts of the 
echo were contaminated by stars. The mean surface brightness 
of the echo is expected to be less vulnerable to the procedure of 
cutting out the stellar patterns. 

We have attempted to correct the measured total flux for this 
efl 'ect. From the radius of the outer rim of the echo measured 
in iTvlenda et al.l ([2005). we calculated the total number of pix- 
els expected inside the echo. Then we defined a covering fac- 
tor, /c, as a ratio of the number of pixels that detected the echo 
(i.e. counted in column 4 of Table [TJ, to the total number ex- 
pected from the size of the echo. This factor is generally be- 
tween 0.7-1.0. Assuming that /c is < 1.0 because of the field 
stars' patterns and that the mean surface brightness is statisti- 
cally the same over the echo image, the corrected echo flux can 
be calculated as the measured flux (listed in column 5) devided 
by /c. The resulting value should, however, be treated as an up- 
per limit to the echo flux. The missing pixels do not result from 
the field stars only, but also from the fact that there were re- 
gions of the echo, not contaminated by stars, where the echo 
was physically too faint to be measured. Therefore, as a finally 
corrected flux, we took a mean value from the flux estimated in 
the above way and the measured flux (column 5). For all except 
the observations from December 17, both fluxes were taken with 
the same weight when averaging. The echo images obtained on 
December 17 show empty regions significantly more extended 
than on earlier dates, so we took the measured flux with a dou- 
ble weight when calculating the finally corrected flux for this 
date. The echo fluxes corrected with the above proc edure and 
converted to the ST magnitudes (ISirianni et al. I I2OOI are given 
in column (7) of Table [1] The corrections to the echo total flux 
were largest for the F435W filter, but even in this case they do 
not exceed 10% (0.10 mag.). 

As noted above, the covering factor, /c, is generally < 1 .0, but 
it varies from image to image. It is usually lowest for the F435W 
filter and shows a tendency to decrease with the time of obser- 
vations. As a result the surface brightness listed in column (6) 
of Table [T] is not entirely comparable between different filters 
and diff'erent dates. This may create uncertainties and ambigui- 
ties when comparing the results of the echo simulations to the 
observational measurements. Therefore we also derived another 
mean surface brightness, which was calculated from the bright- 
est pixels in each image, counted until they fill /c = 0.5. The 
resulting values, named 5b 05 below, are Usted in column (8) of 
Table [U 

Table [1] lists uncertainties of the obtained values in addi- 
tion to the results of measurements. The errors of the total flux 
and the surface brightness (columns 5, 6, and 8) are expressed 
in percent of the values, while those of STmag (column 7) are 
given in magnitudes. When evaluating the errors we took into 
account count statistics, uncertainities due to subtracting the sky 
background, uncertainties of the echo dimension, and calibra- 
tion errors. The first source of errors, which we evaluated from 



the standard deviation of the mean sky background, was negli- 
gible compared to the other ones and resulted in errors ^0. 1 % 
in the flux and surface brightness, and $0,001 mag in STmag. 
As noted in Sect. |2] before measuring the echo, we clipped the 
pixels with signal below 1.5cr of the sky variations. This value 
is of course somewhat arbitrary. To estimate the uncertainties 
caused by this procedure we repeated measurments by allowing 
for Icr variations of the clipping level. The surface brightness, 
Sb,, was most vulnerable to this uncertainty, which was 4-7% 
in the F606W and F814W filters, and 8-16% in the F435W fil- 
ter The total flux was uncertain to 0.3-1.3% in the F606W and 
F814W filters and 2-8% in the F435W filter. 5bo.5 is insensi- 
tive to the clipping level. Errors in determining the echo radius 
derived in iTvlenda et al.l ( l2005h propagate through the correc- 
tion procedure (described above) when deriving STmag, result- 
ing in a ~0.05 mag uncertainty. They also affect S'bo.s through 
the evaluation of the total number of pixels inside the echo and 
result in 2.5-3.5% uncertainty. The deriviation of S-r does not 
involve the echo radius. To evaluate uncertainties of the image 
flux calibration, we selected ten field stars located outside the 
light echoes and measured their corresponding counts in all im- 
ages. Calculating mean values and standard deviations, we found 
that the calibration is uncertain to 7%, 5%, and 4% in the F435W, 
F606W, and F814W filters, respectively. Propagation of all these 
errors through the measuring procedure results in the final errors 
listed in Table [T] 

4. The model 

Basic assumptions and formulae o f the hght echo geometry 
and its modelling can be found in ITvlendal (|2004{) (see also 
ISparks et al.l |2008 [). Briefly, we assume a single-scattering ap- 
proximation and that the echo structure is described by the 
paraboloid equation, i.e. 

r^Z + ct, (1) 

where r is a radial distance of the scattering point from the source 
of radiation, z is a projection of r on the line of sight of the 
source, and f is a time delay between the moment of observa- 
tions of the scattered radiation and that of the radiation directly 
recorded from the source. 

The geometry assumed in our light-echo simulations is 
shown in Fig. [1] The echoing medium surrounding the flaring 
object is uniform and has a semi-infinite geometry. Its boundary 
is in a form of a plane perpendicular to the line of sight, situated 
in front of the flaring object at a certain distance, zo, from the 
object. The medium extends well beyond the object so that the 
light echo paraboloid never reaches its possible boundary behind 
the object during the time span covered in our simulations. It was 
shown in Tvlenda (2004) that V838 Mon is situated in a cavity in 
the surrounding dusty medium. We assume, for simplicity, that 
this dust-free hole is spherical in respect to the object and has a 
radius ro. 

In the above approximation, the echoing medium is 
parametrized, apart from zo and ro, with the optical thickness 
of the medium along the line of sight of the central object, tq. 
The thickness devided by zo - ro gives the extinction coeffi- 
cient, gext = 2abs + fisca, where 2abs and gsca are the absorption 
and scattering coefficients, respectively. The scattering on dust 
grains is usually anisotropic, so in order to model it properly, we 
adopted a standard anisotropy phase function, i.e. 



(1 +g2_2g cos 0)3/2' 
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Table 1. Measurements of V838 Mon's light echo on images obtained with HST/ACS. 



date 


obs.mode 


total 


pixels 


flux 


surf.brightness 


STmag 


surf.brightness(/(. = 0.5) 


2002 




count s"' 


number 




Wm - ' 


arc sec - 


corrected 


Wm^^ytim"' arcsec"^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


(8) 


F435W 


Apr.30 


polarimetric 


9.59E4 


4.20E5 


2.95E-13 ±7.3% 


2.81E-16d 


-AQ.1% 


12.72 ±0.07 


4.14E-16±7.5% 


May 20 


polarimetric 


7.55E4 


4.80E5 


2.33E-13±8.6% 


1.94E-16d 


:13.0% 


12.92 ±0.08 


2.65E-16±7.4% 


Sept.2 


standard 


3.58E4 


9.44E5 


l.lOE-13+7.8% 


4.67E-17d 


:12.0% 


13.69 ±0.08 


6.36E-17±7.5% 


Oct.28 


standard 


2.90E4 


1.21E6 


8.92E-14 ±7.3% 


2.96E-17d 


:10.3% 


13.94 ±0.07 


4.16E-17±7.5% 


Dec. 17 


standard 


2.48E4 


1.27E6 


7.64E-14±11.2% 


2.41E-17d 


:17.5% 


14.10±0.11 


3.10E-17±7.4% 



F606W 



May 20 polarimetric 1.02E5 6.28E5 7.90E-13 ±5.0% 5.82E-16 ±7.5% 11.66 ±0.05 9.46E-16 ±5.8% 

Sept.2 polarimetric 5.19E4 1.16E6 4.01E-13 ±5.0% 1.41E-16 ±7.3% 12.39 ±0.06 2.25E-16 ±5.8% 

Oct.28 standard 3.94E4 1.36E6 3.05E-13 ±5.0% 8.93E-17 ±7.0% 12.68 ±0.05 1.44E-16 ±5.8% 

Dec. 17 polarimetric 3.52E4 1.45E6 2.72E-13 ±5.0% 7.48E-17 ±7.2% 12.77 ±0.05 1.09E~16 ±5.7% 



F814W 



May 20 standard 1.29E5 5.07E5 8.98E-13 ±4.4% 7.08E~16 ±8.9% 11.49 ±0.05 1.16E-15 ±5.6% 

Sept.2 standard 7.04E4 1.07E6 4.89E-13 ±4.3% 1.83E-16 ±6.8% 12.14 ±0.05 2.87E-16 ±5.3% 

Oct.28 standard 5.96E4 1.36E6 4.14E-13 ±4.3% 1.22E-16 ±5.8% 12.34 ±0.04 1.97E-16 ±5.2% 

Dec. 17 standard 5.40E4 1.51E6 3.75E-13 ±4.5% 9.95E-17 ±8.4% 12.43 ±0.05 1.50E-16 ±5.1% 
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Fig. 1. Cross section, y = 0, of the geometrical model adopted in 
our light-echo modelling. The source of the light is at x = 0, 
y - 0, z = (black point in the figure). The observer is 
at jc = 0, y = 0, z = oo. jc and z are in units of ct. Dust 
is uniformely distributed in the space satisfying z < Zo and 
r = -sjx^ + z^ > ro. Dashed curves: light-echo paraboloids, 
Eq. ([U, for two epochs differing by a factor of 2. If the source 
produced a flash of light, which is observed between these two 
epochs, dust scattering the light at present (ct = 0) towards the 
observer is situated between these two paraboloids. 



where is a scattering angle and g = (cos 0} is an anisotropy 
factor. Generally, the anisotro py factor is d etermined by prop- 
erties of dust grains (see e.g. lDrain3ll985h . In our modelling, 
(cos 6*) is treated as a free parameter to be determined from fitting 
the results of simulations to the observational measurements. 
Extinction of the light travelling through the echoing medium 
is taken into account in the modelling. 



The observed light curve of V838 Mon was obtained in the 
Johnson-Cousins filters. The HST images of the light echo were 
obtained with the HST filters, however, which differ from the 
Johnson-Cousins system. For the purpose of our modelling we 
converted the photometric measurements of V838 Mon in the 
B, V, Rr, and Ic, taken from the same references as in Tylenda! 
UOm, into the fight curves in F435W, F606W, and F814W us- 
ing transformation formulae of lSirianni et al. I (l2005l) . These light 
curves were then used to model the evolution of the light echo 
images. 

As discussed in Sect. |2] the observed echo image was con- 
taminated by field stars, in particular by the diffraction pattern of 
V838 Mon itself. The latter had a typical dimension of ~3 arc- 
sec in radius (note that the diffraction cross was much more 
extended), and was removed in the reduction precedure of all 
images before measuring the echo. This effect can be easily 
taken into account in the model simulations and we removed the 
central circular region of the above radius from calculating the 
model echo flux and surface brightness; 

Foflowing Sparks et al. (2008) andl Afsar & Bondl(l2007h . we 
assume that V838 Mon is at a distance of 6.1 kpc. For a fixed 
distance, the expansion of the echo outer rim is uniquely deter- 
mined by Zo - T he observed ex pansion of the light echo outer rim, 
as measured in'Tvlenda et al.' (^2005) and compared to the results 
of our modelling in Fig. |2l results in zo - 2.1 pc. Following 
Tylenda (2004) and the adopted distance, we assume ro =0.1 pc. 
Note that the latter value is not crucial for the results of our anal- 
ysis. The albedo of dust grains, i .e. Qsca/Qext, was taken from 
the interstellar extinction curve of lWeingartner & Draind (1200 ih 
withRv = 3.1. 

5. Results of model fitting to observations 

The results of our modelling with the fixed values of the distance, 
Zo, '"0, and albedo, as described in Sect.|4] depend on two param- 
eters, i.e. optical thickness of the echoing matter along the line of 
sight in front of V838 Mon, To, and the scattering anisotropy fac- 
tor, (cos 9). This is illustrated in Fig. [3] Note that each of the two 
parameters affects the model results in a different way. As can be 
seen from the upper panel of Fig. [3] to does not change the shape 
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Fig. 2. Observ ed expansion of the outer rim of the light echo, 
as measured in lTvlenda et alj (l2005h (symbols), compared to our 
modelling (lines) assuming a distance to V838 Mon of 6.1 kpc. 
The curves were obtained taking zq -2.1 pc (middle curve) + 
10% (outer curves). The outer rim in the simulations was defined 
as the radius at which the echo surface brightness drops to 10% 
of its mean value. The light curve of V838 Mon, used in the 
simulations, is the same as that used to model the light echo in 
the F435 W filter, see upper panel of Fig. |4] 



of the echo flux evolution but it moves the curve vertically in the 
diagram. On the other hand, (cos 6) determines the rate of fad- 
ing of the light echo with time (see the bottom panel of Fig. |3j. 
Therefore fitting the model results to the observed evolution of 
the light echo in different filters can constrain the values of these 
two parameters and their dependence on the wavelength. 

Figures|4]and|5]compare the results of our models of the light 
echo (curves) fitted to the observational measurements (sym- 
bols with error bars) collected in Table [T] The fitting was made 
with the method. Figure |4] compares the total brightness of 
the light echo in the ST magnitude scale (data from column 7 
in Table [T]l. Here, we also plotted the observed light curves of 
V838 Mon (dashed curves) that were used to model the echo. 
Figure |5]presents the observed and model evolution of the mean 
surface brightness of the echo. Open symbols and full curves re- 
fer to the measurments and modelling of 5 b (data from column 
6 in TablellJ. Filled symbols and dashed curves shows the same 
but for 5bo.5 (data from column 8 in Tablellli. The parameters of 
all models displayed in Figs. |4] and |5] are given in Table |2] The 
errors of the parameters were obtained form the 90% confidence 
level of the;^'^ fitting. 
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Fig. 3. Evolution of the model echo total flux (in the magnitude 
scale) as a function of the optical thickness, tq, along the line 
of sight of the central source (upper panel) and the anisotropy 
factor, (cos 6), (bottom panel). Dashed curve in both panels: the 
light curve of the central source (the same as in the upper panel 
of Fig.|4]i. Full curves in the upper panel: the light echo flux evo- 
lution with (cos 6} fixed at 0.6 but for tq varying from 0.05 (bot- 
tom curve), through 0.1, 0.2, 0.3, and 0.5 (uppermost curve). Full 
curves in the bottom panel: the light echo flux evolution with 
To fixed at 0.2 but for (cos 6} varying from 0.0 (bottom curve), 
through 0.2, 0.4, 0.6, and 0.8 (uppermost curve) 



6. Analysis and discussion 

As can be seen from Table |2] the optical thickness of the dusty 
matter in front of V838 Mon, tq, as determined from modelling 
the surface brightness, ^b, (middle rows in the table) is system- 
atically higher than the values derived from fitting the total echo 
flux (upper rows in the table). The difference is within the error 
bars, however Modelling of the surface brightness of the bright 
regions covering 50% of the echo surface, 5 b 0.5, (bottom rows in 
the table) gives significantly higher values of tq than in the two 
previous cases. This is not surprising, as brighter regions mean 



more effective scattering and thus a higher extinction coefficient. 
The values listed in the upper and middle rows of Table |2] can 
be considered as typical for the echoing medium of V838 Mon. 
Those in the bottom rows, when compared to the other ones, il- 
lustrate how inhomogeneous the medium is. 

The results in Table [2] allow us to study properties of dust 
grains in the echoing medium. Figures|6]and|7]show values of the 
scattering coefficient and the anisotropy factor, respectively, as 
functions of wav elength and compare th e m to t he analogous val- 
ues obtained by IWeingartner & Draind (12001 1) from modelling 
the interstellar extinction curve with Ry - 3.1 and 5.5. 
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Fig. 4. Evolution of the light-echo total-flux in the ST magnitude 
scale as observed in the F435W (upper panel), F606W (mid- 
dle panel), and F814W (bottom panel) filters. Dashed curve: ob- 
served light curve of V838 Mon. Asterisks: observed light-echo 
magnitudes (column 7 in Table[T]i. Full curve: best fit of the mod- 
elled evolution to the observational points. Parameters of the fit 
can be found in Table|2]upper rows. 
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Fig. 5. Evolution of the light-echo surface-brightness in the 
F435W (upper panel), F606W (middle panel), and F814W (bot- 
tom panel) filters. Open symbols and full curve: observed light- 
echo surface -brightness (column 6 in Table [T]) and its best- 
fitting model. Filled symbols and dashed curve: the same but for 
Sbq.5 (column 8 in Table [TJ. Parameters of the fit can be found 
in Table 12] middle and bottom rows. 
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Table 2. Values of tq and (cos 9) derived from the model fitting 
to the observations. 



F435W 



F606W 



F814W 



total flux (STmag) 










To 0.204 ± .016 
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.003 


0.087 ± 


.003 


(cos G) 0.67 ± .06 


0.62 ± 


.02 


0.59 ± 


.03 


surface brightness 5 b 










To 0.2 18 + .028 


0.159 ± 


.009 


0.087 ± 


.006 


(cos(9> 0.64 ±.10 


0.57 ± 


.06 


0.54 ± 


.06 


surface brightness 5b os 










To 0.306 + .020 


0.249 ± 


.010 


0.137 ± 


.004 


(cos e) 0.64 + .05 


0.60 ± 


.04 


0.56 ± 


.03 
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Fig. 6. Dependence of the scattering coefficient, Qsca, on the 
wavelength as derived from modelling of the echo. The val- 
ues are normalized to Qsca at the effective wavelength of the 
F814W filter Asterisks: results of the models fitting the total 
flux, STmag, (upper rows in Table|2]l. Open circles: results of the 
models fitting the mean surface brightness, 5b, (middle rows in 
Table|2]i. Filled circles: results of the models fitting the mean sur- 
face brightness derived from the brightest pixels filling - 0.5, 
5'bo.5, (bottom rows in TableO. The horizontal error bars repre- 
sent the widths of the photometric bands. Open and filled sym- 
bols are shifted by +0.05 yum in the abscissa axis for clarity. 
Curves: relations expected from modelling the interstellar ex- 
tinction curve of Weingartner & Draine (2001) for Ry - 3.1 
(dotted) and Ry - 5.5 (dashed). 



As can be seen from Fig. |6l our modelling of different ob- 
servational quantities (total flux, surface brightness) results in 
a similar wavelength dependence of the scattering coefficient. 
We can also conclude that this relation is not significantly differ- 
ent from that expected from ffie standard composition and size 
distribution of dust grains used to model the interstellar extinc- 
tion curve of Weingartner & Draine (2001) with Ry - 3.1. This 
suggests that the dusty matter responsible for the observed echo 
is not significantly diff'erent from the typical interstellar dusty 
medium. 

The scattering anisotropy factor, (cos 9), as displayed in 
Fig. I2I is more sensitive to the observational uncertainties than 
the scattering coefficient. There are systematic differences in the 
derived (cos 9} values, depending on which quantity was anal- 
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Fig, 7. Dependence of the anisotropy factor, (cos 9), on the wave- 
length, as derived from our modelling of the echo. The symbols 
and curves have the same meaning as in Fig.|6] 



ysed, but the differences are within the error bars. The value of 
(cos 9} is determined from the rate of the echo fading with time. 
This rate, if derived from the measurements of all pixels detect- 
ing emission from the echo, is likely to be distorted because the 
portion of the echo surface covered by these pixels is different 
for different dates of the observations. Therefore we expect that 
the results of fitting the mean surface brightness, Sso.5, (bottom 
row in Tableland filled circles in Fig.|7|i are more reliable than 
the other two. 

Fig. Q shows that our determinations of (cos 9} gave values 
systematically greater than those expected from the interestel- 
lar extinction curve with Ry - 3.1. They are better reproduced 
by the extinction curve with Ry - 5.5, but even in this case the 
points mostly lie above the curve. Greater (cos 9} is expected for 
bigger dust grains, but this usually implies that the scattering co- 
efficient is less dependent on the wavelength, which would be in 
conflict with the points in Fig.|6] However, the model extinction 
curve and (cos 9} depend not only on the grain size, but also on 
the chemical composition and structure of the grains. It is thus 
possible that playing with different parameters of the dust grains 
would allow one to better reproduce the dust properties we ob- 
tained from the analysis of the light echo of V838 Mon. This is 
out of the scope of the present study, however. 

We can estimate from Table |2] (upper and middle rows) that 
the additional extinction of V838 Mon due to local dust respon- 
sible for the light echo is ^ 1.2 ro(F606W) ^ 0.18 (Ee-y ^ 
0.05). Of course, this is the most likely value derived from the 
mean properties of the local dusty medium seen in the echo. 
This medium is inhomogeneous, however, as is shown by the un- 
even distribution of the surface brightness over the echo images. 
Therefore the real extinction along the line of sight of V838 Mon 
can be different from the above estimate. Nevertheless, the re- 
sults obtained from the brighter half part of the echo surface, i.e. 
from fitting Sbo.s (bottom rows in Table |2]i, indicate that it is 
rather improbable that Ay > 0.3. With the observed extinction 
curve (Fig.|6]l this would also imply Eg-y ^ 0.10 from the local 
dust, which would be too large compared to the observed dif- 
ference in Es-y between V 838 Mon and the thre e B-type main- 
sequence stars observed bv lAfsar & Bondl (l2007i) (see Sect.[TJ. 
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Using the sta ndard interstellar-m edium relation, Ay/Nu ^ 
5.3 X 10"^^ cm^ (iBohlin etaT]ll978h . where A^h is the column 
density of hydrogen (atomic plus molecular), and taking Ay - 
0.18, we obtain a surface density of the echoing matter in front 
of V838 Mon to be ~ 8 x lO""* g cm^^. In 2004 the angular ra - 
dius of the light echo reached ~1 arcmin jTvlenda et al.ll2005l) . 
With a distance of 6.1 kpc this allows us to estimate the mass 
of the diffuse matter of the light echo to be ~35 Mq. This is a 
lower limit to the mass of the circumstellar matter of V838 Mon, 
since, first, in later epochs the echo expanded to even larger 
radii, second, the estimate is based on the optical thickness in 
front of V838 Mon only, while the matter almost certainly ex- 
tends beyond the object. The above estim a te can be compared to 
90-100 Mo obtained bv iKaminski et all (1201 Ih for the molec- 
ular complex seen on the CO rotational lines in the vicinity 
of V 838 Mon and a hundred M© estimated bv iBaneriee et al.l 
(|20"06) for the mass of the echoing matter in the infrared. All 
these estimates clearly show that the diffuse matter in the vicin- 
ity of V828 Mon, partly seen in the light echo, could not have 
resulted from a past mass loss activity of V838 Mon, as advo- 
cated, e.g. in Bond et al. (2003) and Bond (2007). It is much 
more probable that we see remnants of an interstellar complex, 
from which V838 Mon, its companion, and perhaps the other 
members of the observed cluster were formed. 

It is worth noting that Kerve lla et al.l (1201 2l ) recently found 
that the mass of the dusty matter echoing pulsating radiation 
from the Cepheid RS Pup is ~290 Mq. It seems therefore that, 
in general, a stellar light echo to be well resolved in imaging 
observations requires a flaring star to be associated with a dusty 
medium of interstellar origin rather than formed from a mass 
loss activity of the star 

On the basis of our results, we can quite safely conclude that 
the optical thickness of the dusty matter seen in the light echo 
is too low to account for the observed dimming of the B-type 
compani on by ~1.3 mag compared to the other B-type stars ob- 
served bv lAfsar & Bondl(l2007i) in the vicinity of V838 Mon (see 
Sect. [1]). Other explanations of this discrepancy have to be con- 
sidered, e.g. mis-classification of the B-type companion, an evo- 
lutionary status of the companion different from that of the three 
B-type stars of the cluster, or a mistaken distance to V838 Mon 
and its companion. 

The classificati on of the companion as a B3 main-sequence 
star com es from Munari et al.l (l2002 b'. see also Muna ri et al.l 
(I2007hl)1. Thi s clas sification was essentially confirmed by 
lAfsar & Bondl (|2007|) . but the B-ty pe component in their spec- 
trum, similarly as in the spectrum of^ Kamiriski et al.1 (12009 ). was 
then significantly contaminated by V838 Mon. The observed V 
brightness of the companio n would be i n accord with those of 
the three B-type stars of Afsar & Bondl (12007) if it were of a 
~B7 V spectral type. The intrinsic colours of su ch a star would 
be (B - V)q = -0.13 and (U - B)q = -0.43 (ISchmidt-Kalerf 
Il982h . From Goranskij's' compilation of the photometric data 
we derived V = 16.19 ± 0.03, B - V = 0.55 + 0.05, and 
U - B - -0.15 + 0.07 for the B-type companion. The observed 
B -V colour can be reconciled with that of the above standard 
if Eb-v - 0.69. But then the reddened U - B colour of the stan- 
dard star is -1-0.08, i.e. it is significantly redder than the observed 
one. Usin g the photometric me asurements of the B-type com- 
panion of iMunari et al ] (I2007bh . we obtain Eb-v - 0.82 and 
U - B - -1-0.17 for the reddened standard, which is to be com- 
pared to the observed one of -0.06. We thus conclude that B7 V 
cannot be reconciled with the observed colours of the B-type 
companion. Note that the best fit of a standard star to the colours 
derived from the Goranskij data can be obtained taking the spec- 



tral type of B2.5 V and Eb-v - 0.80, which is no t significantly 
different from the B3 V spectral type obtained by iMunari et al.l 
(20023). 

The second possibility is that the B-type companion of 
V838 Mon is intrinsically fainter because it is in a some- 
what different evolut ionary stage than the three B-type stars of 
lAfsar & Bondl (l2007h FI All existing observational data indicate 
that all four stars are on the main sequence (MS) or not far from 
it. The only reasonable situation to be discussed in this context 
is that the B-type companion of V838 Mon has a similar effec- 
tive temperature butj^s less luminous (by ~1.3 mag.) than star 9 
(B3 V) of lAfsar & Bondl (i2007l) because it is less massive but 
is at its hottest stage, namely near the zero age main sequence 
(ZAMS), while star 9 is more evolved on the MS. Taking the 
calibrations of the absolute magnitude versus spectral type and 
luminosity class from Schmidt-Kaler (ll982|), one easily con- 
cludes that to account for a ~ 1 .3 magnitude difference between 
two B-type stars of a similar spectral type with one of them be- 
ing on ZAMS, the second star has to be of luminosity class IV, 
i.e. slightly off MS. This would probably not be in a significant 
conflict with the spect roscopically de t ermine d class V, given the 
quaUty of the data of lAfsar & Bondl (|2007|) . However, the dif- 
ferences in magnitude between star 9, star 8 (B4 V), and star 7 
(B6 V) imply that the latter two stars have also to be of luminos- 
ity class IV. To have a situation in which three stars of different 
masses are about to leave or just left MS at the same time, one 
has to conclude that the three stars have different ages, e.g. star 7 
would have to be about twice as old as star 9. This disproves the 
idea that these three stars form a cluster. Futhermore, a case in 
which three B-type stars formed at different times happen to lie 
in a close vicinity in the sky and mimic an isochronic sequence 
of a cluster in the colour-magnitude diagram is extremely im- 
probable. So is the whole scenario discussed in this paragraph. 

If two MS stars of the same spectral type suffer from a 
similar extinction but are significantly different in the observed 
brightness then the only reasonable conclusion is that the stars 
are at differe n t distan ces. In this way we return to the study of 
IMunari et al.l (l2007al) . where the authors concluded, primarily 
on the basis of the spectroscopic distance of the B-type com- 
panion, that V838 Mon is at a distance of ~10 kpc. In this case 
V838 Mon and its companion would have nothing to do with the 
cluster of Afsar & Bond (2007). This possibility creates a prob- 
lem with the observ ed Unear pola rization of the light echo, how- 
ever, as analysed in ISparks et al.l (i200 8). To reconcile the latter 
with the distance of 10 kpc, one has to postpone that maximum 
polarization occurs at a scattering angle of ~60° and not at 90°, 
as assumed in Sparks et al. (2 008J. Theoreti cally this cannot be 
excluded (e.g. Chapter 4.2 in lKruegelll2008l) but requires rather 
unusual dust grains, e.g. metallic particles. 

We performed light echo simulations, similar to those de- 
scribed in Sect.jS] but assuming a distance of 10 kpc. They gave 
results not significantly different from those obtained in Sect. |5] 
except that (cos 9) was systematically higher by ~20% than the 
values in Table |2] 

We can conclude that if the distance of V838 Mon were 
~10 kpc, dust in the echoing medium would have to be quite 
peculiar, i.e. giving an extinction curve close to the standard one 
but scattering strongly in forward directions and giving maxi- 
mum polarization at scattering angles significantly lower than 
90°. With the present state of our knowledge of interstellar dust 



' As sugested by V. Goranskij during the STScI workshop 
"Intermediate-Luminosity Red Transients", June 2011. 
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grains, we consider the above result as an argument against a 
distance as large as ~10kpc. 

In summary, we cannot present any strightforward sce- 
nario that could explain the ~1.3 magnitude difference between 
the B-type compani on of V838 Mon and the B3 V star of 
lAfsar& Bondl Qm% and which would be in accord with the 
observed properties of dust in the echoing matter of V838 Mon. 
Perhaps this could be a dust cloud giving a grey extinction, ly- 
ing well in front of V838 Mon (not to be seen in the light echo), 
and not intervening wi t h the lines of sight of the three B-type 
stars of I Afsar & Bondl (l2007h . A detailed spectroscopic study, 
particularly of the interstellar line profiles (Nal, KI), as well 
as multi-colour photometric measurements of a large sample of 
stars in the field of V838 Mon would perhaps enable one to be 
conclusive in this subject. 
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